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Abstract 

Recently, Yakubovich [Opuscula Math. 26 (2006) 161-172] and Passian et al. [J. Math. 
Anal. Appl. doi:10.1016/j.jmaa.2009.06.067] have presented alternative proofs of an orthogo- 
. nality relation obeyed by the Macdonald functions of imaginary order. In this note, we show 



that the validity of that relation may be also proved in a simpler way by applying a technique 
occasionally used in mathematical physics to normalize scattering wave functions to the Dirac 
delta distribution. 
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CN , 1 Introduction 
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I In recent papers, Yakubovich [1] and Passian et al. [2] proved the following orthogonahty relation 

. ■ for the Macdonald functions of imaginary order: 

o> 

dx '-^ '-^ ^ ————S{iy - ly') {iy,i^'>0). (1.1) 

J X sinh(7ri^j 

^"-^ ' 

In the first of these works, a proof rehed on advanced techniques of the theory of distributions. 
■ An approach adopted in the second paper was a two-step one. At first, two heuristic arguments 

JH \ making the relation (jl.ip plausible were presented. One of these arguments was based on an integral 

relation between the Macdonald functions of imaginary order and Mehler's conical functions, for 
which a counterpart orthogonality relation had been known for a long time. The second argument 
given in support of the validity of Eq. (jl.lD exploited the fact that the Laplace transform of Kii,{x) 
is a known elementary function. Subsequently, a sophisticated proof of the relation (jl.ip , different 
from the one in Ref. [1], was presented. 

It is the purpose of this note to present still another proof of the relation ()1.1|) . The approach 
we adopt here is known in mathematical physics, where it is occasionally used for normalization 
of scattering states to the Dirac delta distribution. In some sense, it is akin to the standard 
method used to prove weighted orthogonality relations for eigenfunctions of regular Sturm-Liouville 
problems (cf, e.g., Ref. [3, Sec. 7.1]). It has several advantages. First, it is elementary compared to 
the methods used in Refs. [1,2]. Second, it is constructive: one derives the orthogonality relation. 
Finally, it is general and may be used to obtain counterpart orthogonality relations not only for 
Kiy(x), but also for other special functions (e.g., Ref. [4]). 
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2 A summary of relevant properties of the Macdonald func- 
tions of imaginary order 

Before we proceed to the merit, in this short section we shall summarize these properties of the 
Macdonald functions of imaginary order which will be exploited later in Sec. [3l The formulas 
presented below have been excerpted from the collection of Magnus et al. [5] . 

The function Kii,{x) is a particular solution to the modified Bessel differential equation 



2d^F{x) AF{x) 



+ {v"^ - x^)F{x) = 0. 



dx^ dx 

Other particular solutions to Eq. (|2.ip are the modified Bessel functions of the first kind 
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k=0 
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The relationship between the three functions is 

Kiy[x) = — — — ^ , 

2i smh(7ri^j 

from which it follows immediately that 

Ki^[x) = K^i^{x). 

For large positive values of x, the function Ki^{x) has the asymptotic representation 



/^e-^[l + 0(.-)], 



while for x — > +0 from Eqs. (|2.3p and (|2.2p and with the aid of the known relationship 



one deduces that 



KUx) 



V sinh(7ri/) 



-i^ In — + argr(ii/) 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



O X sin 



-i^ln- + argr(2 + ii^) 
(i/eR). (2.7) 



3 Derivation of the orthogonaUty relation for the Macdon- 
ald functions of imaginary order 

To derive the orthogonality relation for the Macdonald functions of imaginary order, we proceed as 
follows. If i4rii/(x), with 1/ G M, is substituted for F{x) into Eq. (|2.ip . this results in the differential 
identity 

d / di^i^(a;)\ /t^^ \ , , „ 

^ X ! + x\K,,{x)=^. (3.1) 



da; \ da; ^ 

The counterpart identity for the function K\yi(x)^ with v' G M, is 



d (^ dK,Ax) \ ^ ( :__^[ j^.^,(^)^o_ (3 2) 



dx V da; 
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Next, we premultiply the first of the above equations by Kii,i{x), the second one by Ki^{x), subtract 
and integrate the result over x from x = ^>Otoa; = oo. After obvious rearrangements, this gives 



ax \ ax J ax \ ax 



The integral on the right-hand side of Eq. (|3.3|) is easily evaluated by parts; one obtains 

oo 

Ki^{x)Ki^'{x) 



1/2 -j.'2) / dx 



X \ Kiy[x) — Kii,,{x) 



dx 



(3.3) 



(3.4) 



x=5 



By virtue of Eq. (12. 5|) . the expression in the bracket on the right-hand side of the above relation 
vanishes in the upper limit. Hence, we obtain 



^ KUx)K,.'{x) /"^^^^ de ""^^^ de 
1 = - ^ 73 77^ 



(3.5) 



and consequently 



da; 



Kii,(x)Ki^,{x) 



lim f- 



(3.6) 



Use of the asymptotic representation (|2.7|) and of elementary trigonometric identities transforms 
Eq. jHSl) into 



da; 



K--,i,{x)Ki^>{x) 



2^J vv' sinh(7ri/) sinh(7ri'') 
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— (y — v') In ^ + argr(izy) — argr(iz/') 



sm 



-(i/ + i/') In I + argr(ij/) + argr(ii/') 



(3.7) 



To evaluate the limit on the right-hand side of Eq. (|3.7p . we observe that if f{ri) is a real analytic 
function of 77 G R, such that /(O) = (which implies that lim,,^o /('y)/'? is finite), then in the 
distributional sense it holds that 

a+f{ri)/n a oo 

hm ""[°'^ + ^^"^] = hm / da e'- = ^ lim / da e'- = / da e'-. (3.8) 

a^oo TTTy ZTT a^oo J ZTT a^oo J 27r J 



In the expression at the extreme right of the above chain of equalities one immediately recognizes 
the well-known Fourier representation of the Dirac delta distribution S{r]) (cf, e.g., [6, Sec. 4.5]), 
so that, provided the function /(yy) satisfies the above constraints, one has 



lim ^"^N + /(^)U ^(^). 

a-+oo nrj 



(3.9) 
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As for +0 it holds that — ln(^/2) oo, with the help of the above relationship Eq. (|3.7p 
becomes 

oo 

dx ^^'^^^^^"^'^^^ = -^=^L===[6i. - + 6i. + .')]. (3.10) 
X z^/t^t^ sinh(7rz/) sinh(7rzy'j 

Exploiting in Eq. (I3.10|) the following basic property of the delta distribution [6, Sec. 4.4]: 

9{v')S{v-r/)=g{v)S{v-r/), (3.11) 

one eventually arrives at the sought orthogonality relation 

oo 

/ dx = - — — — r (5 -iy ) + S{L' + i')]. 3.12 

J X 21/ sinh(7rzy) 



If we impose the constraint i', v' > 0, then v + v' > Q and consequently in the distributional sense 
we have 

5{u + u')^Q {u,u'>0). (3.13) 

It is then evident that under the above restriction the orthogonality relation (|3.12p turns into the 
one in Eq. p.ip . 
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